Observation of standing kink waves in solar spicules 
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Abstract We analyze the time series of Call H-line 
obtained from Hinode /SOT on the solar limb. The 
time-distance analysis shows that the axis of spicule 
undergos quasi-periodic transverse displacement at dif- 
ferent heights from the photosphere. The mean period 
of transverse displacement is ~180 s and the mean am- 
plitude is 1 arc sec. Then, we solve the dispersion re- 
lation of magnetic tube waves and plot the dispersion 
curves with upward steady flows. The theoretical anal- 
ysis shows that the observed oscillation may correspond 
to the fundamental harmonic of standing kink waves. 
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1 Introduction 

Observation of oscillations in solar spicules may be 
used as an indirect evidence of energy transport from 
the photosphere towards the corona. Transverse mo- 
tion of spicule axis can be observed by both, spec- 
troscopic and imaging observations. The periodic 
Doppler shift of spectral lines have been observed from 
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ground based coronagraphs (INikolskv fc Sazanovlll987 ; 
Kukhianidze et allliooS IZaaarashvili et al.ll2007l ). But 
Doppler shift oscillations with period of ~5 min also 
have been observed on the SOlar and H eliospheric Ob- 
servatory (SOHO) bv IXia et aD (|2QQ5h . Direct peri- 



odic displacement of spicule axes have been found by 
imaging obse rvations on Optical Sola r Telescope ( SOT) 
on_J^ dej|De Pontieu et al.l 120071 : iKim et all 12008 : 
He et al. 2009[ ). The torsional Alfven waves were re- 
ported recently in the context of a flux tube connecting 
the photosphere and the chr omosphere as per iodic vari- 
ation of spectral line width ( Jess et al. 2009[ ). 

The observed transverse oscillations of spicul e axes 
were interpreted by kink (INikolskv fe Sazanov 1987 : 



Kuk hianidze et a l. 2006; Zaqarash vili et al. l2QQ7 tlKim et alj 
20081 ) and Alfven ([De Pontieu et al J l2QQ7h waves. All 



spicule osc illations events are summarize d in a recent 
review by Zaqarashvili fc Erdelvi ( 2009[ ). They sug- 
gested that the observed oscillation periods can be for- 
mally divided in two groups: those with shorter peri- 
ods (<2 min) and those with longer periods (^2 min) 
( Zaqarashvili fc Erdelvil l2009h . The most frequently 
observed oscillations lie in the period ranges of 3-7 min 
and 50-110 s. 

Spicule seismology, which means the determina- 
tion of spicule properties from ob served oscillations 
and w as originally suggested by IZaqarashvili et al 



(2007)) , has been significantly developed during last 
years (|Aiabshirizadeh et~aL 20091 : Verth et al. 2011 : 
Tavabi et al]l201ll ). 

Spicules are almost 100 times denser than sur- 



rounding coronal plasma ([Beckers! 1 19681 ), therefore 
they can be considered as cool magnetic tubes em- 
bedded in hot coronal plasma. Wave propagation 
in a static m agnetic cylinder was st udied by Edwin 
and Roberts ([Edwin fc Roberts! fl983l ) . They derived 
general dispersion relation of all possible wave modes 
in magnetic tubes. Then the linear and non-linear 
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MHD waves propagation in a cylindrical magnetic flux 
tube with axial steady flo ws have been also studied 



( Terra- Homem et all 120031 ) . They show that steady 



flows change the treatments of propagating waves be- 
cause of induced Doppler shifts. 

In the present work, we study the observed oscil- 
lations in the solar spicules through the data obtained 
from Hinode. Then we model the oscillations as magne- 
tohydrodynamic (MHD) waves in magnetic flux tubes. 



in respect with the hypothetic line (white line on the 
images), which is drown on the same place during whole 
time series. The clear quasi-periodic transverse motion 
of the spicule axis is seen on the figure. 

We use the time slice diagrams in order to study the 
quasi-periodic motion of the spicule axis in detailed. 
Figure [2] shows the time slice diagrams performed at 6 
different heights from the limb. Each cut is obtained 



2 Observations and image processing 

We used a time series of Call H-line (396.86 nm) 
obtained on 25 October 2008 during 03:20 to 03:25 
UT by the Solar O ptical Telescope onboard Hinode 
( Tsuneta et~aDl2008l ). Note, that the Call H-line ob- 
servations of the same day have b een used recently t o 
study multi-component spicules bv lTavabi et al.l ([201lf h 
but they used another time interval of these data. The 
spatial resolution reaches 0.2 arc sec (150 km) and the 
pixel size is 0.1082 arc sec (~80 km) in the Call H-line. 
The time series has a cadence of 10 seconds with an 
exposure time of 0.5 seconds. The position of X-center 
and y-center of slot are, respectively, arc sec and 948 
arc sec, while, X-FOV and Y-FOV are 112 arc sec and 
56 arc sec respectively. 

We used t he "fgprep," "fgrigidalign" and "madmax" 
algorithms ( Koutchmv fc Koutchmv 1989| ) to reduce 



the images spikes and jitter effect and to align time 
series and to enhance the finest structures, respectively. 

Figure [T] shows 9 selected images of the time series, 
which consists of 27 consecutive images. We study the 
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Fig. 1 9 images of time sequence in Call H-line (left to 
right and top to bottom). White diamonds indicate the 
spicule locations in time series. 

transverse motion of selected spicule (black narrow line) 




Fig. 2 Time slice images of time series shown in Fig. [T] 
The x-axis corresponds to the heliocentric coordinates mea- 
sured from the white diamonds of Fig. [T] 

by averaging over 9 pixels along the spicule axis, which 
corresponds to 1 arc sec around each height. 

The dark regions represent the spicule at particu- 
lar height. We clearly see that the spicule axis un- 
dergoes the transverse oscillation at each height with 
nearly same periodicity. The oscillation period at each 
heights is estimated as ~180 s. The oscillation am- 
plitude is nearly 1 arc sec, but slightly changes with 
height. Figure [3] shows the height variation of the os- 
cillation amplitude. The amplitude increases almost 
linearly with height and its slope is 0.06. The oscilla- 
tions of spicules axis are closed to the standing pattern 
as we do not see any upward or downward propagation. 

T he density is al most homogeneous along the spicule 
axis ([Becker si 1 19681 ). therefore the density scale height 
should be much longer that the spicule length. Then the 
variation of oscillation amplitude with height can not 
be due to the decreased density. Therefore, we assume 
that the height dependence of amplitude is due to the 
standing oscillations. Then we may have that the oscil- 
lation amplitude is proportional to sin(/cz), where k is 
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Fig. 3 The height variations of displacement amplitude of 
the spicule axis oscillation. 



the oscillation wave number. This expression can be ap- 
proximated in the long wavelength limit as ~kz. Then, 
the oscillation wavelength is A ~ 1/0.06 arcsec ~ 11500 
km. This allows the estimation of the kink speed as 
v k = X/T ~ 11500/180kms -1 - 65 kms -1 . The esti- 
mated period and kink wave speed are in good agree- 
ment with the period and speed of fundamental har- 
monic of kink waves. 

The oscillation amplitude and phase speed allows to 
estimate the energy flux, F, storied in the oscillations 



F = -pv v k , 



(1) 



where p, v and Vk are density, wave amplitude and 
phase speed respe ctively. The d ensity in spicules is 
3 x 10 -10 kgm~ 3 (|Beckerd[l968l ). The wave velocity 
can be determined as v = £o/^o, where £o is the axis 
displacement estimated as 1 arcsec and To is the os- 
cillation period estimated as ~180 s. We take the es- 
timated phase speed as Vk = 65 kms -1 . With these 
parameters the energy flux is estimated as F = 1500 



J m 



Total coronal energy losses in the quiet Sun 



is 3000 Jm s , therefore the energy flux storied in 
the oscillation is of the same order as the energy losses, 
but probably is not enough to fully compensate them. 



3 Theoretical modeling 

We use ideal MHD equations to model the propaga- 
tion of waves in spicules. The density seem s almost 
unifo rm along almost whole length of spicules ([Beckers 
19681 ), therefore we ignore the effects of stratification 
due to gravity. The density inside spicules is almost 



100 times larger than outside, therefore we model them 
as magnetic flux tubes embedded in coronal magnetized 
plasma. We consider a uniform vertical magnetic tube 
of radius a with uniform magnetic field Bqz (B e z) inside 
(outside) the tube. The kinetic gas pressure and den- 
sity inside (outside) of the tube are po (p e ) and po (p e ), 
respectively. We also consider uniform steady flows in- 
side (outside) the tube with velocity Uoz (U e z). The 
pressure balance condition at the tube boundary im- 
plies that 
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The densities inside and outside the tube are related as 
([Edwin fc Robertslll983t ) 
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(3) 



where c VlPo/po (c e = \/lVe/pe) and v A o = 
Bo / ^/ppo (vAe — B e /y/pp e ) are the sound and Alfven 
speeds inside (outside) the tube, respectively. Here 7 is 
the ratio of specific heats and p is the magnetic perme- 
ability. 

Fourier transform of linearized MHD equations with 
the assumption that all the perturbed quantities are 
ocexp[i(— Lit + nO + kz)] and the contin uity of the per- 
turbed interface ([Chandrasekharl Il96lh and the total 
pressure across the cylinder boundary r = a yield the 
dispersion relation ([Terra- Homem et al.l 120031 ) 
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(4) 



where I n and K n are modified Bessel functions of or- 
der n. Here m a and Q a (a = 0, e) are given by the 
expressions 



m2 _(^-W"L-^) 



(4 + «L)(* 2 4a-n|) 



and 



Q a = uo — kU a , 



where 



c a v Aa 



(5) 



(6) 



(7) 



is the tube speed. This dispersion relation describes 
both, surface (rag > 0) and body waves (mQ < 0). 

We solved the dispersion relation for the values of 
c s o = 10 kms -1 , c se = 262 kms -1 , vao = 40 kms -1 , 



^Ae = 210 kms -1 and po/Pe = 0.02. With these data 
we have Ck = 42.5 kms -1 , cto = 9.7 kms -1 and 
c Xe = 100.6 kms -1 . Figure H] shows the dispersion 
diagram of kink waves under above conditions in the 




Fig. 4 The phase speeds of the kink modes (the curve 
with green color is fundamental kink mode) under spicule 
conditions. 

absence of steady flows. All speeds in the plots are nor- 
malized to ^aq. As VA e > vaq then the surface waves 



are absent ([Edwin fc Roberts! 1 19831 ) . Figure [5] shows 



the same dispersion diagram but in the case of steady 




Fig. 5 The same as in Fig. [H but for upward flow of 25 
kms -1 . 

flow Uo = 25 kms -1 . It is seen that phase speeds and 
consequently wave frequencies are shifted due to the 
steady flows. 

Let us assume that the full length of spicule is L. 
Then the wavelength of the fundamental standing mode 
is A = 4L and thus the wave number of the fundamental 
mode is k = 2tt/X = 7r/2L. The second and the third 
harmonics have the wave numbers tt/L and 2tt/L , re- 
spectively. The mean length of classical spicules varies 
from 5000 to 9000 km in H a ([Zaaarashvili fc Erdelvil 
l2Q09h . Then the periods of standing kink modes can 



be estimated from the dispersion diagram. For the the 
spicule length of 3500-9000 km we obtain the period of 
fundamental, second and third harmonics in the ranges 
of 170-320 s, 85-160 s and 40-80 s, respectively. The 
observed oscillation period of ~180 s may correspond 
to the fundamental harmonic of standing kink waves. 



4 Discussion and conclusion 

We performed the analysis of Ca II H-line time series 
at the solar limb obtained from Hinode/ 'SOT in order 
to uncover the oscillations in the solar spicules. We 
concentrate on particular spicule and found that its 
axis undergos quasi-periodic transverse displacement 
about a hypothetic line. The period of the trans- 
verse displacement is ~ 180 s and the mean ampli- 
tude is ~ 1 arc sec. The sa me periodicity was found 
in Doppler shift oscillation b y lZaqarashvili et al.l (j2007l ) 
and lDe Pontieu et al.l (j2007l h so the periodicity is prob- 
ably common for spicules. 

We model spicules as dense plasma jets along mag- 
netic flux tube being injected from the chromosphere 
upwards into the hot corona. Solving the wave disper- 
sion relation in the presence of steady flows we conclude 
that the steady flows change the characteristics of wave 
propagation in a cylindrical magnetic flux tube. The 
calculated periods of fundamental, second and third 
harmonics of standing kink modes with an upward flow 
of 25 kms -1 are in the ranges of 170-320 s, 85-160 
s and 40-80 s, respectively, for the spicule length of 
3500-9000 km. 

Therefore, the observed quasi-periodic displacement 
of spicule axis can be caused due to fundamental stand- 
ing mode of kink waves. The energy flux storied in the 
oscillation is estimated as 150 J m -2 s -1 , which is of the 
order of coronal energy losses in quiet Sun regions. 
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